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ABSTRACT 

The Boston University-Five College Radio Astronomy Observatory (BU- 
FCRAO) Galactic Ring Survey (GRS) of 13 CO J = 1 -> emission covers Galac- 
tic longitudes 18° < t < 55.7° and Galactic latitudes \b\ < 1°. Using the SE- 
QUOIA array on the FCRAO 14m telescope, the GRS fully sampled the 13 CO 
Galactic emission (46" angular resolution on a 22" grid) and achieved a spectral 
resolution of 0.21 km s -1 . Because the GRS uses 13 CO, an optically thin tracer, 
rather than 12 CO, an optically thick tracer, the GRS allows a much better deter- 
mination of column density and also a cleaner separation of velocity components 
along a line of sight. With this homogeneous, fully-sampled survey of 13 CO 
emission, we have identified 829 molecular clouds and 6124 clumps throughout 
the inner Galaxy using the CLUMPFIND algorithm. Here we present details of 
the catalog and a preliminary analysis of the properties of the molecular clouds 
and their clumps. Moreover, we compare clouds inside and outside of the 5 kpc 
ring and find that clouds within the ring typically have warmer temperatures, 
higher column densities, larger areas, and more clumps compared to clouds lo- 
cated outside the ring. This is expected if these clouds are actively forming stars. 
This catalog provides a useful tool for the study of molecular clouds and their 
embedded young stellar objects. 

Subject headings: catalogs-molecular data-ISM:clouds-ISM: molecules-Galaxy: 
kinematics and dynamics 
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Introduction 



The Boston University-Five College Radio Astronomy Observatory (BU-FCRAO) Galac- 
tic Ring Survey (G RS) is a survey of G alactic 13 CO J = 1 — ► emission toward the first 
Galactic quadrant (jjackson et al. 2006 ). One of the primary goals of the GRS was the 
study the Milky Way's so-called 5kpc molecular ring. Discovered as a peak in the radial 



Galactic CO distribution midway between the sun and the Galactic center ( iBurtonl 11976 
Scoville fc Solomon!ll975l ). the 5 kpc molecula r ring contains ~ 70% of the molecular mate 



rial within the solar circle (M ~ 2 x 10 9 M Q ; ICombeslll99ll ). It is the prominent feature in 
molecular maps of the Galaxy and is the location of the majority of Galactic star formation 
jBurtonl Il976f iRobinson et al.lll984l : Idemens et"aJll988h . 



Because the 5 kpc ring covers such a large region of sky, a large, high angular resolution, 
sensitive survey is required to reveal the molecular clouds and dense clumps where star 
formation occurs. In order to be able to study their properties, the BU-FCRAO GRS was 
designed to map the 13 CO molecular line emission over a large portion of the 5 kpc ring 
and, in doing so, answer a number of fundamental questions that remain about the 5 kpc 
ring. For instance, how many molecular clouds does the 5 kpc ring contain? What are their 
properties? Do molecular clouds in the star-forming ring differ from clouds elsewhere in the 
Galaxy? 

The GRS provides a unique database to identify star-forming molecular clouds and 
their dense clumps. The first step in such studies is to isolate, identify, and catalog the 
molecular clouds and their clumps. Combined with distance determinations, this catalog 
can characterize the cloud masses, sizes, line widths, and densities in a range of Galactic 
environments. This paper describes a catalog of molecular clouds and clumps identified 
using 13 CO data from the GRS which provides a new, homogeneous database for the study 
of th e large- and small -scale structure of molecular gas throughout the inner Milky Way 



Simon et al. 2001 ). The complete catalog is available in electronic forrrfl In total, 



we 



(e.g. 

identify 829 clouds and 6124 clumps. A preliminary analysis of their general characteristics 
is also presented. In particular, we compare the properties of those clouds located inside the 
5 kpc molecular ring to those outside the ring. 



see http: / /www. bu.edu/galacticring/ molecular_clouds.html 
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2. The BU-FCRAO Galactic Ring Survey 



The BU-FCRAO GRS Jjackson et all2006h mapped Galactic 13 CO J = 1 -> emission 
using the SEQUOIA multi-pixel array on the FCRAO 14 m telescope. The survey extended 
in Galactic longitude from £ = 18° to 55.7° and in Galactic latitude from \b\ < 1°, for a total 
of 75.4 square degrees. The survey's velocity coverage was —5 to 135 km s -1 for Galactic 
longitudes £ < 40° and —5 to 85kms~ 1 for Galactic longitudes £ > 40°. The typical rms 
sensitivity was cr(T* A ) ~ 0.13 K. The survey comprises a total of 1,993,522 spectra. The data 
are available to the community at www.bu.edu/galacticring or in DVD form by request. 

Unlike most previous surveys of the inner Galaxy, the GRS is fully sampled (46" an- 
gular resolution on a 22" grid), has a high spectral resolution (0.21 kms -1 ), and traces 
13 CO J = 1 — > emission. Compared to 12 CO, the 13 CO molecule is ~ 50 times less abun- 
dant and, thus, has a much lower optical depth. As a result, 13 CO is a much better tracer 
of column density. Moreover, because it suffers less from line blending and self-absorption, 
it can separate clouds along the same line of sight more easily than 12 CO. Because of these 
improvements over previous surveys the GRS can detect many new structures and clou d 
cores previously missed by older 12 CO surveys (e.g.. ISanders et al.lll986l ; iDame et al.ll200ll ). 



For more details of the telescope and instrumental parameters, the observing modes, 
the data reduction processes, and the emission and noise characteristics of the dataset see 



Jackson et al.l (120061 ). 



3. Identifying molecular clouds and clumps 



We use the automated algorithm CLUMPFIND (j Williams et al .1119941 ) to identify molec- 
ular clouds and clumps within the GRS. CLUMPFIND searches through a three-dimensional 
(£, b, v) data cube using iso-brightness surfaces to identify contiguous emission features with- 
out assuming an a priori shape. In this respect, it is more flexible than other cloud identi- 
fication algorithms that artifically decompose the emission features into three-dimensional 
Gaussian profiles (e.g., GAUSSCLUMPS). 

CLUMPFIND begins the search for clouds at the voxel with the peak brightness in 
the data cube. The algorithm steps down from this peak brightness in levels referred to as 
'contour increments'. During the first iteration the algorithm finds all contiguous voxels with 
brightnesses between the peak value and the next level (defined as the peak value minus the 
contour increment). If these contiguous voxels are isolated they are assigned to a new cloud. 
If they are contiguous with a previously identified cloud then they are assigned to that pre- 
existing cloud. The algorithm iterates until it reaches a minimum brightness level. At this 
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lowest level, CLUMPFIND does not search for new clouds, but instead adds all emission at 
that level into pre-existing clouds. This procedure is straightforward for isolated emission 
peaks, but complicated for blended emission feature s. For blended emiss ion, features are 
separated based on a 'friends-of-friends' algorithm (see lWilliams et al.lll994j for more details). 

For each cloud (or clump) that is identified, CLUMPFIND gives the position of the 
peak of the emission for the cloud in (£, b, v), the peak temperature, the FWHM extent of 
the cloud in (I, b, v), an estimate of the radius of the cloud (assuming the total number of 
voxels form a sphere), the total number of voxels, and the sum of the emission within the 
voxels. Moreover, a three-dimensional (£, b, v) data cube is also created which is identical 
to the input data cube in (I, b, v) space. In this output cube, however, the temperature 
scale in each voxel is replaced by an integer which corresponds to the cloud to which it is 
assigned. This output cube is extremely useful for identifying the exact voxels associated 
with each individual cloud and is used extensively in the analysis described below. 



3.1. Clouds 

Giant molecular clouds (GMCs) are generally defined as extended molecular line emis- 
sion features with typical properties as outlined in Table [TJ Molecular clouds do not have 
uniform densities, but are observed to be fragmented on all size scales. We refer to the 
dense regions within GMCs where clusters may form as 'clumps' and the smaller, denser 
sites where the individual star formation occurs as 'cores'. Table [1] compares the typical 
properties of each of these structures. 



3.1.1. Smoothing of the GRS data 

To improve the signal-to-noise and to better identify faint, extended molecular clouds we 
smoothed the GRS data both spatially (to a resolution of 6') and spectrally (to a resolution 
of 0.6 km s _1 ) using Gaussian kernels. Because this smoothing was moderate, it did not cause 
a significant loss of information in £, b, or v. In particular, this smoothing did not cause 
molecular clouds to be artifically blended along the line-of- sight. The smoothing improved 
the rms sensitivity from ~0.13K for the original data to ~ 0.01 K for the smoothed data. 
The smoothed data were then re-sampled on a 3' grid. Figured] shows a comparison between 
the original and smoothed data in both (£, b) and (£, v) space. 
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3.1.2. Modification of the CLUMPFIND algorithm 

The CLUMPFIND algorithm was originally designed to identify the clumpy substruc- 
ture within individual molecular clouds. To run CLUMPFIND, the user must specify two 
input parameters; the contour increment and the lowest brightness level. In the original 
algorithm, the lowest brightness level is hardwired to be a multiple of the contour incre- 
ments which are forced to be evenly spaced, from the peak brightness to the minimum level. 
Because of the large dynamic range of the GRS data, however, the algorithm in this original 
form was unable to identify simultaneously both extended, diffuse clouds as well as compact, 
bright clouds. Moreover, it could not easily separate clouds that appeared distinct, espe- 
cially at the lowest emission levels. For instance, by specifying a large contour increment 
(e.g. 0.7 K) CLUMPFIND would identify the brightest clouds in the GRS as one contiguous 
feature, but missed or merged many of the fainter, more extended clouds. However, with 
a smaller contour increment (e.g. 0.1 K), CLUMPFIND would better identify and separate 
the fainter, extended clouds, but the largest, brightest clouds would be artificially dissected 
and identified as numerous separate clouds. 

Thus, in order to identify molecular clouds within the GRS at all emission levels, we 
modified the CLUMPFIND algorithm. For this purpose we required a large contour incre- 
ment, but a small minimum brightness level. To achieve this, we changed the minimum 
brightness level to be a user input variable independent of the contour increment. In this 
scheme, we retained evenly spaced contour levels between the peak and minimum bright- 
ness levels as per the original algorithm but could additionally search for clouds to a lower 
brightness level. This scheme allowed us to identify both faint, extended clouds and bright, 
compact clouds with the same input parameters. 

3.1.3. Selection of CL UMPFIND input parameters 

With this modified version of CLUMPFIND, we tested a range of contour increments and 
lowest brightness levels to determine which set of parameters would best identify molecular 
clouds. Because of the inherent difficulties in determining where one cloud ends and the 
next begins, different choices of these parameters will certainly result in different output 
catalogs. However, regardless of the exact parameters, the bright, large molecular clouds are 
typically always recovered; the differences in the catalogs is most noticeable in the separation 
of nearby emission features and for the faint, small clouds. 

To test the performance of the algorithm we completed a series of simple tests to compare 
the clouds identified with CLUMPFIND to the original dataset. In essence we replaced the 
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CLUMPFIND clouds with elliptical three-dimensional Gaussian models with the same peak 
temperature, center position, size, orientation, line width and voxels of each cloud identified 
by CLUMPFIND. We then subtracted these Gaussian models from the original dataset and 
examined the residuals. We describe this process below. 

We ran three representative regions within the smoothed GRS (each covering A£ x Ab 
= 2°x2°) through CLUMPFIND, varying both the contour increment and the minimum 
brightness level. The test regions were selected to span a range from crowded and complex 
to sparse and simple emission features. We varied the contour increment from 0.1 (10cr) to 
0.7 K in steps of 0.1 K and the minimum brightness level from 0.05 (5a) to 0.4 K in steps 
of 0.05 K and examined the output from all combinations of these two parameters. 

Comparisons between the Gaussian model and the original dataset reveal that, for 
simple unblended emission features, the Gaussian model matched the data very well: the 
mean residual between the data and the model was ~ 2a of the data. Not surprisingly, the 
mean residuals were smallest when the minimum values of the contour increment and lowest 
minimum brightness level were used (0.10 K and 0.05 K respectively). Put another way, the 
Gaussian model reproduced the data best when a larger number of clouds were found (e.g., 
a residual ~0.023 K for N c i ou d s = 115 compared to a residual ~0.028 K for N c i ou d s = 9). 

In regions of bright, complicated emission features and blended lines, the mean residual 
between the original dataset and the Gaussian model were slightly higher ( residual ~ 0.065 
K). For these regions, the mean residuals were insensitive to the choice of CLUMPFIND 
input parameters. For the range of CLUMPFIND input parameters tested, the residuals 
varied by only ~ 0.01 K (la). Thus, for identifying the brightest molecular clouds in the 
dataset (> 10a) we found the exact choice of parameters is not crucial. 

Figure [2] shows an example of the identification of an isolated molecular cloud in (£, b) 
and (£, v) space. Because of the difficulties of displaying three-dimension data cubes in a 
two-dimensional image, we have selected to show a single velocity channel in the (£, b) image 
(left panels) and a single Galactic latitude plane in the (£, v) image (right panels). The top 
panels show an isolated molecular cloud in the smoothed GRS data (in color scale and as 
contours in this and all subsequent images). The white ellipses mark the approximate extent 
of the cloud as determined by CLUMPFIND; the major and minor axes are equivalent to 
the projected extent of the cloud in each direction. The middle panels of Figure [2] show 
the corresponding Gaussian model, which was generated using the peak temperature, center 
position, size, orientation, and line width output from CLUMPFIND. The lower panels show 
the CLUMPFIND output cube; the color scale in these images represent the voxels that are 
assigned to each cloud for each particular channel or plane. 
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Williams et al.l (119941 ) found that CLUMPFIND most accurately represents the data 
when the contour increment is set to twice the 1 a noise of the data. Because we are 
interested in selecting the largest and brightest molecular clouds in our dataset, we use 
a contour increment of twice the 10a noise (rather than la). Thus, to identify clouds 
we use a contour increment of twice the 10a noise in the smoothed data (0.20 K) and a 
lowest brightness level of 20a (0.20 K). Our comparisons of the Gaussian models with the 
original data indicated that these parameters did not produce extremes in either the number 
of detected clouds or mean residuals. Thus, we consider them suitable for isolating and 
identifying molecular clouds. 

The output from CLUMPFIND was also checked by eye against the original dataset 
to verify that these parameters could simultaneously identify and separate both bright and 
compact molecular clouds in addition to diffuse and extended molecular clouds. This choice 
of input parameters was confirmed to be satisfactory and were adopted to produce the final 
catalog. Thus, using a contour increment of 0.20 K and a lowest brightness level of 0.20 K 
CLUMPFIND, identified 848 molecular clouds in the GRS dataset. 

For clouds to be selected as real features, CLUMPFIND requires each cloud to have a 
minimum number of voxels (we set this number to 16). To be confident we are identifying 
real molecular clouds, we imposed additional size criteria on each cloud before inclusion in 
the final list. We exclude all clouds that have their measured sizes in any axis less than 
or equal to the resolution of the smoothed data in that axis (i.e., A£ or Ab < 6' or A V 
< 0.6 km s" 1 ). We find 18 clouds that meet these criteria. Thus, when we exclude these, 
the final catalog contains 829 molecular clouds. In the sections below we characterize the 
properties of these clouds. 



3.I.4. Determining sizes of the clouds 

Because molecular clouds have complicated morphologies, it is difficult to quantify their 
shapes and measure their sizes. Although CLUMPFIND determines the extent of each 
cloud in Galactic longitude and latitude, this determination fails to account for the cloud's 
orientation. To more accurately describe their shapes, we modeled the two-dimensional 
integrated emission for each cloud using a two-dimensional ellipse fitting routing. Two- 
dimensional integrated intensity maps of each cloud were produced by summing in velocity 
over all the voxels identified by CLUMPFIND for a given cloud. The center position of the 
cloud was then determined from the centroid of this velocity integrated image (with all pixels 



3 written in IDL by D. Fanning. 
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weighted equally regardless of their integrated intensity). The semimajor and semiminor axes 
were determined by fitting an ellipse to this integrated intensity image. The listed position 
angle was measured counter-clockwise from the positive Galactic longitude axis. With this 
convention, a cloud whose major axis lines along the Galactic plane will have a position 
angle of 0° while one lying perpendicular to the Galactic plane will have a position angle of 
90°. This method provides a simple description of the cloud's position, size, and orientation 
with respect to the Galactic plane. 

We have also calculated the projected two-dimensional area for each cloud which was 
determined by multiplying the total number of pixels within the integrated intensity image 
by the area of a pixel in square degrees. 

3.2. Clumps 

Because molecular clouds show structure on all size scales, they typically contain several 
compact, dense clumps. These clumps presumably represent the compact, dense regions 
within the molecular clouds where star and cluster formation may take place. To identify 
these star-forming clumps, we searched for three-dimensional substructure within each of the 
829 molecular clouds in (£, b, v) space. To achieve this we used the full angular and spectral 
resolution GRS data and the original CLUMPFIND algorithm. 

For each molecular cloud, we search for clumps using only the voxels for that cloud 
identified by CLUMPFIND as described in the previous section. We also restricted the size 
of each clump to at least 50 voxels (in total in (I, b, v) space) in order to prevent false 
detections. 

The same technique using three-dimensional Gaussian models (as described in §3.1.31) 
was used to determine the best CLUMPFIND input parameters to identify and separate the 
clumps within each of the clouds. We selected 10 molecular clouds to test the CLUMPFIND 
input parameters for clump identification. These molecular clouds were selected to include 
a range in sizes, emission levels, and substructures. To test the clump identification, we 
varied the contour increment using values of 0.2, 0.26, 0.36, 0.39 K and the lowest brightness 
level from 0.52 K to 1.3 K in steps of 1.3 K and examined, by eye, the output from all 
combinations of these parameters. 

We found a contour increment of 0.26 K (2a) and a minimum brightness level of 1.30 K 
(10cr) was best at isolating and identifying the individual clumps within the clouds. Using 
these parameters, we identify, in total, 6135 clumps within these clouds. Again, for the final 
catalog, we exclude all clumps that have sizes smaller than the resolution in any axis (i.e. 
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A£ or Ab < 0?01 or A V < 0.2kms 1 ). Thus, the final clump catalog contains 6124 entries. 



3.3. Examples of clouds and clumps 

Figure [3] shows the integrated intensity images for two clouds identified in the GRS. The 
(£, b) images were produced by integrating the emission over their velocity range, such that 
the integrated intensity I = J T mb dv. Because it is very difficult to represent the asymmetric, 
three-dimensional output of CLUMPFIND in a two-dimensional image, we have included on 
these images white ellipses that indicate the approximate area of the region identified by 
CLUMPFIND for each clump. To indicate the extent of the clouds, we use two-dimensional 
ellipses whose major and minor axes are equal to the projected extent of the clump in each 
direction. The identification clumps is straightforward for GRSMC G018. 14+00.39, but 
more complicated for GRSMC G043.34-00.36 (see Fig. EJ. 

We find the number of clumps within these 829 molecular clouds range from 1 to 111, 
with a typical molecular cloud containing ~ 7 clumps. Because some clouds are too faint, we 
detect no clumps within them above the clumpfind threshold. However, most of the clouds 
(~ 96%) contain at least one clump that has a peak T m & > 10 times the average temperature 
of the cloud. These are the brightest clumps in the clouds where star-formation will likely 
occur. 



4. Description of the catalog 



Table [2] gives a sample of the catalog entries for each molecular cloud. Table [3] gives a 
sample of the catalog entries for the clumps. The complete catalogs are available in electronic 



format at http://www.bu.edu/galacticring/molecular_clouds.html 



4.1. Clouds 

For each molecular cloud identified, we list parameters output from the CLUMPFIND 
algorithm in addition to several derived quantities. The columns of Table [2] are as follows: 
(1) the molecular cloud name, designated as GRSMC (for GRS Molecular Cloud) followed 
by the Galactic longitude and latitude (t and b) coordinates of the peak of the emission in 
degrees, e.g. GRSMC G053. 59+00.04; (2), (3) and (4) the Galactic coordinates {£ and b) and 
velocity (Vlsr) of the peak voxel of the emission; (5) the velocity FWHM defined as 2.35 
times the velocity dispersion (A V); (6) main beam temperature of the peak voxel (T m b); 
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(7) and (8) the Galactic coordinates {t and b) of the centroid of the integrated intensity 
image; (9) and (10) the semimajor and semiminor axes (a and b), determined from the 
elliptical Gaussian fit to the integrated intensity image; (11) position angle of the fitted 
ellipse, measured counter-clockwise from the positive longitude axis (PA); (12) the projected 
two-dimensional area (A); (13) the average main beam temperature of the cloud, defined as 
the average main beam temperature of all the voxels (T„„); (14) the peak integrated intensity 
(Ipeafc)j (15) the total 13 CO integrated intensity (It taJj)j (16) the peak H 2 column density 
N(H 2 ); and (17) a flag denoting if the cloud lies on one of the survey boundaries ('X', 'Y' or 
'V indicate that the cloud lies on the boundary in Galactic longitude, Galactic latitude or 
velocity, respectively) 

Because molecular clouds have irregular shapes and complex morphologies, the parame- 
ters listed in columns (7)— (11) are given simply to estimate the clouds' approximate position, 
size, and orientation. 

4.2. Clumps 

The clumps listed here represent the compact, dense substructures within the clouds. 
Because we are interested in selecting the brightest, densest regions within each of the clouds, 
not all clouds have entries in this list; for the faintest clouds, we often detected no clumps. To 
describe the shape of the clumps we simply use the extent in Galactic longitude and latitude 
output directly from CLUMPFIND (rather than the elliptical fitting that was performed for 
the clouds). 

The columns of Tableware as follows: (1) the molecular cloud name (from Table [2]); 
(2) the clump number, e.g. cl, c2, c3, etc; (3), (4) and (5) the Galactic coordinates (£ and 
b) and velocity (Vlsr) of the peak voxel of the emission; (6), (7) and (8) the FWHM extent 
in longitude, latitude and velocity (A£, Ab, and AV); (9) the main beam temperature of 
the peak voxel (T mb ); (10) the projected two-dimensional area (A); (11) the peak integrated 
intensity (l pea k)'i (12) the total 13 CO integrated intensity (Jtotai)] (13) the peak H 2 column 
density (N(H 2 )); and (14) a flag denoting if the clump lies on a boundary. 



4 where hotai = f J T m & dv dfl 
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5. Ensemble properties 

5.1. Galactic distribution 

Figure H] shows the number distribution in Galactic latitude (left) and longitude (right) 
of the 829 molecular clouds identified in the GRS. Included in this figure for comparison is 
the integrated intensity image for the GRS, that is the 13 CO emission integrated over all 
velocities. The location of the clouds in (£, v) space is shown in Figure El overlaid on the 
GRS (£, v) image, which is the 13 CO emission averaged over all Galactic latitudes. The 
clouds identified span the complete range in Galactic longitude, latitude, and velocity as 
covered by the GRS. 

Because the Milky Way rotation is well approximated by an axisymmetric rotation 
curve, the line of sight velocity directly relates to a Galactocentric radius. Thus, the (£, v) 
positions can uniquely determine Galatocentric radii. To determine Galactocentric radii, 



we assume circular motions and the rotation curve of I Clemens! (119851 ) with (Ro,6o) = (8.5 
kpc, 220kms~ 1 jE Because we are sampling clouds primarily in the flat part of the rotation 
curve, the derived kinematic distances are insensitive to the exact choice of rotation curve. 
Figure M shows the number distribution of clouds as a function of Galactocentric radii. The 
peaks in this distribution correspond to the known spiral features within the Galaxy. 

In the inner Galaxy, a single Galactocentric radius corresponds to two kinematic dis- 
tances, one at the near distance and one at the far distance. Accurate distances cannot 
be determined until this ambiguity is resolved. A comparison of GRS data with Hi self- 
absorption can resolve this ambiguity and will be presented in a future paper (J. Duval et 
al. in prep.) 



5.2. Properties of the clouds and clumps 

Figure [7| shows the number distributions of the clouds (diagonally hashed histograms) 
and clumps (horizontally hashed histograms) as a function of peak T m b, line width, semima- 
jor axis, semiminor axis, position angle, and peak H2 column density. 

We find that the clouds have a mean peak 13 CO T m & of ~ 1.6K, AV of ~ S.Gkms^ 1 , 
semimajor axes of ~ 0?41, and semiminor axes of ~ 0?23. The fact that the semimajor axes 
are typically twice the semiminor axes suggests that most clouds are elongated. In fact, the 



5 Note that the Clemens rotation curve explicitly calls for a small velocity correction factor due to a 
measured mis-calibration of the local standard of rest. 
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mean value of the ratio of semimajor to semiminor axis is 2.0. On the other hand, we find 
that the clumps are typically brighter (peak T m b of ~ 5.2 K), have lower line widths (AV 
~ 1.4kms _1 ), are smaller (semimajor axes of 0?06), and are rounder (the semimajor and 
semiminor axes are comparable). 

The position angles measured for the clouds range from —90 to 90° with a peak value 
of V.2 (Fig. [7]). The position angles were measured counter-clockwise from the positive 
Galactic longitude axis, so a peak value in the distribution of 1?2 implies most clouds are 
orientated along t he Galactic plane. This result is in agreement with previous studies using 
the GRS dataset (IKoda et al.ll2006l ). This result may arise from an observational bias due 
to the extended areal coverage of the GRS along the Galactic plane (~ 38° were covered 
in Galactic longitude compared to the 2° covered in Galactic latitude). The fact that few 
clouds were detected toward the edges of the coverage in Galactic latitude, however, suggests 
that our results are not severely effected by this lack of latitude coverage. 

The peak H2 column density, N(H2), was calculated by assuming optically thin 13 CO 
emission with an excitation temperature of 10K and using standard conversion factors (see 
Simon et al.ll200ll ) via the expression 



N(H 2 ) = 4.92 x 10 M I peak 



(cm 



where I pea k is the peak integrated intensity (Kkms -1 ). We find that the clouds have a mean 
peak N(H2) of ~ 6.3 x 10 21 cm~ 2 while the clumps have slightly lower peak column densities 
of ~ 2.4 x 10 21 cm -2 . This is opposite to what is expected if the clumps are the densest 
parts of clouds where star-formation will likely take place. However, this results is easily 
explained by the fact that the clump distribution includes all clumps within a cloud, makes 
the mean of the N(H 2 ) distribution lower for the clumps compared to the clouds. 

Table H] summarizes these parameters for the clouds and clumps and lists the minimum, 
maximum, mean, median, standard deviation, and slope of a power-law fit to the distri- 
butions. We find that the clouds identified have smaller average brightness temperatures, 
larger line widths and sizes, and higher peak column densities compared to the clumps. 



5.3. Determining the excitation temperatures and opacities 



To derive the excitation temperatures (T ex ) and opacities (r) of the clouds within our 
catalog, we make use of the 12 CO J = 1 — ► University of Massachusetts- Stony Brook sur- 
vey (UMSB; ISanders et al.lll986l ). The survey region covered the same region of the Galactic 



13 



plane as the GRS, but mapped the molecular clouds via their 12 CO emission. Direct com- 
parisons between the measured 12 CO and 13 CO temperatures are facilitated by the fact 
that the two surveys were obtained with the same telescope (the FCR AO 14 m). However , 
because the UMSB survey was under-sampled ( 45" beam on a 3 / g rid; ISanders et al.lll986l ) 
compared to the GRS (46" beam on a 22" grid; I Jackson et al.ll2006l ). the combination of the 
two dataset can only provide rough average values for the entire cloud. 

Because the 12 CO J = 1 — > emission is typically optically thick, it can be used to 
estimate the kinetic temperature of the gas via the expression 



T k = 5.532 



in 1 + 



5.532 



T 12 + 0.837 



where T 12 is the 12 CO brightness temperature measured from the UMSB survey. The kinetic 
temperature can then be used to calculate the opacities for the optically thin 13 CO J = 1 — > 
transition using 



kT, 



mb 



hv 



hv 

•kT k _ I 



1 bg 



where is 2.7 K and T m & is the 13 CO brightness temperature measured within the cloud 
from the GRS. To calculate T 12 and T 13 we calculate the mean T m {, from all the correspond- 
ing voxels associated with each cloud as defined by CLUMPFIND within the UMSB and 
GRS datasets respectively. The excitation temperature of the gas is then calculated using 



T, 



mb 



1 - e-r) 



by 



Figure [8] shows the histograms of T ex and r for the clouds. We find that the clouds 
typically have T ex of ~ 9 K and a 13 CO r of 0.13 (see Table |4] for a summary). Most clouds 
have low opacities, which confirms the fact that the 13 CO is optically thin and, thus, a good 
tracer of the column density and mass. 



5.4. Properties of clouds inside and outside the 5 kpc ring 



With a large sample of molecular clouds we can now investigate their properties in a 
range of Galactic environments. For instance, does the temperature, line width, and column 
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density of the molecular clouds within active star-forming regions differ from those in more 
quiescent regions? Given that the majority of Galactic star formation occurs within the 
5 kpc molecular ring, one might expect to see differences in the properties of the molecular 
clouds within the ring, compared to clouds that lie outside the ring. 

To test this idea we have separated the clouds identified with the GRS into those that lie 
inside the 5 kpc molecular ring from those that lie outside the ring. We have defined clouds 
that are associated with the peak in the Galactocentric radius (Rg) distribution (Fig. [6]) at 
~ 4.5 kpc as those associated with the ring (i.e. all clouds with 4 kpc < Rg < 5 kpc). Thus, 
the number of clouds within the ring, Nj n , is 206 (25%), while the number of clouds outside 
the ring, N out , is 623 (75%). 

Figure M shows the number of clouds in and out of the ring (solid and open histograms 
respectively) plotted against their values of T ex , r, line width, peak N(H 2 ), area, and number 
of clumps. Included in the top panels of each of these plots is the fraction of clouds within the 
ring in each of the bins. These plots show that clouds within the ring typically have warmer 
temperatures, higher column densities, larger areas, and more clumps compared to clouds 
located outside the ring. The mean opacities appear comparable between clouds inside and 
outside the ring. We also find clouds within the ring have similar mean values of the line 
width compared to those clouds outside the ring. 

In addition to the general cloud properties, we have also performed a K-S test for each of 
the properties to determine if the distribution in the samples inside and outside of the 5 kpc 
ring are derived from the same parent distribution. Table H] lists the results of the K-S tests 
and shows that for all properties, with the exception of position angles, the distributions 
are not derived from the same parent distribution. Thus, there are clear, and significant 
differences in the distributions of the derived properties for the clouds inside and outside the 
5 kpc molecular ring. 

Warmer temperatures and higher column densities are exactly what is expected for 
active star-forming regions. All the clouds have line widths much greater than the thermal 
line width for gas at their derived T^ (for T^= 10 K, the thermal line width for CO is 
0.13 km s -1 ) suggesting that they are dominated by turbulent motions. 

The fact that we see more clumps within the clouds associated with the 5 kpc molecular 
ring also strengthens the idea that these clouds are forming stars. This suggests that clouds 
within the ring are highly fragmented and have many dense, warm regions where the star 
formation can occur, in contrast to the more quiescent, smoother clouds outside the ring. 

While these results suggest that clouds within the ring have different properties to those 
clouds found outside the ring, we need to consider that clouds within the ring suffer more 
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from line blending and their identification as separate features or merged clouds is more 
dependent on the values input into CLUMPFIND. Moreover, the effects of distance also 
need to be considered; we can more easily separate nearby clouds as opposed to those at the 
far side of the Galaxy. Nevertheless, our analysis suggests that the clouds within the ring 
are significantly different from those outside the ring. 

6. Summary 

Using the CLUMPFIND algorithm we have identified a large sample of clouds and 
clumps within the 13 CO J = 1 — > BU-FCRAO Galactic Ring Survey. In total, we identi- 
fied 829 clouds and 6124 clumps. We find the cloud properties are comparable to typical 
Giant Molecular Clouds, while the clumps are similar to the regions where star and cluster 
formation occur. Moreover, it appears that clouds lying within the 5 kpc ring typically have 
warmer temperatures, higher column densities, larger areas, higher densities, higher masses, 
and more clumps compared to clouds located external to the ring. This difference supports 
the idea that star formation occurs within the ring. We note, however, that there are inher- 
ent difficulties in determining cloud properties without the proper consideration of distance 
biases. 

This catalog provides an invaluable tool for studies of molecular clouds. For instance, 
establishing reliable kinematic distances to the GRS molecular clouds and clumps, and to 
their embedded young stellar objects and clusters, is essential in order to determine their 
masses, sizes, distributions, and luminosities. Thus, combining the catalog of GRS molecular 
clouds with IR Galactic plane surveys from IRAS, MSX, 2MASS, and Spitzer, one can 
obtain luminosities, masses, and sizes of stars and clusters embedded within their natal 
molecular material. One can then address the questions: What is is the spatial distribution, 
luminosity function, and initial mass function of the young stars forming both inside and 
outside the ring? How does the star-formation process within the ring differ from that at 
other locations in the Galaxy, such as near the Sun, in nearby spiral arms, and in the outer 
Galaxy? Moreover, the internal structure of molecular clouds, which traces the influence of 
turbulence in the interstellar medium, can also be studied in a wide range of star-forming 
environments. With a large sample of clouds we can also obtain their clump mass spectra 
and study its relation to the stellar initial mass function. These are all important and 
outstanding questions relating to Galactic star-formation. 

This publication makes use of molecular line data from the Boston University-FCRAO 
Galactic Ring Survey (GRS). The GRS is a joint project of Boston University and Five 
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College Radio Astronomy Observatory, funded by the National Science Foundation under 
grants AST-9800334, AST-0098562, AST-0100793, AST-0228993, & AST-0507657. 
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Fig. 1. — An example of the comparison between the full angular and spectral resolution 
GRS data (in color scale) and the smoothed data (as contours) for GRSMC G019. 39— 00.01. 
The top panel shows a single (£, b) channel while the lower panel shows a single plane in 
(£, v) space. The smoothed data were generated by smoothing both spatially (6' on a 3' grid) 
and spectrally (0.6 km s -1 ) the original data. The smoothed data have an rms sensitivity of 
(j(T* A ) ~ 0.01 K. The contour levels are 0.2 to 3.0 K in steps of 0.3 K on the (£, b) image 
and 0.2 to 2.0 in steps of 0.4 K in the (£, v) image. The white ellipses mark the approximate 
extent of the cloud in (£, b) and (£, v) space. 
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Fig. 2. — An example of the identification of an isolated molecular cloud in (£, b) and (£, 
v) space. The left panels show single velocity channel in the {£, b) image, while the right 
panels show a single Galactic latitude plane in the (£, v) image. The top panels show the 
smoothed GRS data in color scale. On all images the contours are the smoothed data, while 
the white ellipses mark the approximate extent of the cloud as determined by CLUMPFIND; 
the major and minor axes are equivalent to the projected extent of the cloud in each direction. 
The middle panels show the corresponding Gaussian model, which was generated using the 
peak temperature, center position, size, orientation, line width, and voxels output from 
CLUMPFIND. The lower panels show the CLUMPFIND output cube; the color scale in 
these images represent the voxels that are assigned to each cloud for each particular channel 
or plane. 
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Fig. 3.— GRSMC G018. 14+00.39 (left) and GRSMC G043.34-00.36 (right) integrated 
intensity images. Because it is very difficult to represent the asymmetric, three-dimensional 
output of CLUMPFIND in a two-dimensional image, we have include on these images white 
ellipses that indicate the approximate extent of each clump, with the major and minor axes 
equivalent to the projected extent of the clump in each direction. The identification of 
clumps is straightforward for GRSMC G018. 14+00. 39, but more complicated for GRSMC 
G043.34-00.36. 
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Galactic Latitude (degrees) Galactic Longitude (degrees) 

Fig. 4. — The Galactic distribution of the 829 molecular clouds identified in the GRS. The 
top panel shows the integrated intensity image for the GRS (square root color scale, from 
to 40 Kkms -1 ). The lower panels show the number distribution of molecular clouds in 
Galactic latitude (left) and longitude (right). 
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Fig. 5. — The (£, v) diagram for the GRS. The color scale represents the 13 CO emission aver- 
aged over Galactic latitude. The ellipses mark the approximate extent in Galactic longitude 
and velocity for each of the 829 clouds identified. 



-22 - 



80 - 




3 4 5 6 7 8 
Galactocentric radii (kpc) 



Fig. 6. — Galactocentric radial distribution of the molecular clouds identified in the GRS. 
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Fig. 7. — Normalized number distributions for the measured parameters T m b, line width, 
semimajor and semiminor axis, position angle, and peak N(H 2 ) for the clouds (diagonally 
hashed histogram) and clumps (horizontally hashed histogram). Position angles were only 
calculated for the clouds. See Table 0] for a summary of these distributions. 
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Fig. 



Number distributions of the excitation temperature, T eX} (left) and the CO 



opacities, r, (right) for the clouds. These were calculated by combining the 13 CO GRS 
emission with 12 CO J = 1 — » emi ssion from the University of Massachusetts-Stony Brook 
survey (UMSB; ISanders et al.lll986l ). The molecular clouds have mean T ex of ~ 9 K and r 
of 0.13. 
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Fig. 9. — Histograms of the parameters for clouds located in the 5 kpc molecular ring (solid 
histogram) and outside the ring (open histograms). The top panel in each plot shows the 
fraction of the number of clouds in the ring for each bin. We find that the clouds within the 
ring typically have warmer temperatures, higher column densities, larger areas, and more 
clumps compared to clouds located outside the ring. This is expected if these clouds are 
actively forming stars. 
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Table 1: Proper ties of Giant Molecular clouds, clumps, and cores ([Goldsmith! Il987 



Cernicharol Il99ll ) . 



Properties 


GMC 


Clump 


Core 


Size (pc) 


20- 


-60 


3- 


20 


o., c 


)-3 


Density (cm -3 ) 


100- 


-300 


10 3 


-10 4 


10 4 - 


-10 6 


Mass (M Q ) 


10 4 - 


-10 6 


10 3 


-10 4 


10- 


10 3 


Linewidth (kms -1 ) 


6- 


15 


4- 


12 


1 


-3 


Temperature (K) 


7- 


15 


15 


-40 


30- 


100 



Table 2: Properties of the molecular clouds identified in the GRS. 



Cloud 


Peak 


Vls.r 


AV 


^ mb 


Centroid 


a 


b 


PA 


A 




Ipeafc 


^-total 


N(H 2 ) 


Flag 


GRSMC 


I 


b 








I 


b 






















(°) 


(°) 


(km s _1 ) 


(km s _1 ) 


(K) 


(°) 


(°) 


(°) 


(°) 


(°) 


(deg 2 ) 


(K) 


(Kkms" 1 ) 


(Kkms -1 deg 2 ) 


(xlO 22 cm- 2 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


G053.59+00.04 


53.59 


0.04 


23.7 


1.99 


5.75 


53.69 


0.01 


0.52 


0.19 


16 


0.27 


1.90 


36.5 


0.95 


1.8 


- 


G029.89-00.06 


29.89 


-0.06 


100.7 


5.09 


5.46 


29.99 


-0.17 


0.44 


0.35 


-46 


0.46 


2.49 


87.6 


3.12 


4.3 


- 


G049.49-00.41 


49.49 


-0.41 


56.9 


9.77 


5.25 


49.57 


-0.39 


0.38 


0.23 


-21 


0.18 


2.18 


193.2 


1.64 


9.5 


Y 


G018. 89-00.51 


18.89 


-0.51 


65.8 


2.80 


5.17 


18.80 


-0.56 


0.31 


0.31 


-47 


0.29 


2.31 


57.5 


1.51 


2.8 


- 


G030.49-00.36 


30.49 


-0.36 


12.3 


4.56 


4.98 


30.66 


-0.39 


0.47 


0.24 


-28 


0.22 


1.70 


17.5 


0.50 


0.9 


- 


G035. 14-00.76 


35.14 


-0.76 


35.2 


5.00 


4.92 


35.22 


-0.78 


0.31 


0.22 


-11 


0.20 


2.36 


59.3 


1.90 


2.9 


Y 


G034.24+00.14 


34.24 


0.14 


57.8 


5.98 


4.81 


34.19 


0.05 


0.50 


0.41 


-3 


0.52 


1.46 


80.5 


3.15 


4.0 


- 


G019. 94-00.81 


19.94 


-0.81 


42.9 


2.81 


4.58 


19.97 


-0.80 


0.42 


0.18 


-7 


0.22 


1.76 


32.1 


1.00 


1.6 


Y 


G023.44-00.21 


23.44 


-0.21 


101.1 


5.75 


4.40 


23.36 


0.02 


0.64 


0.48 


29 


0.63 


1.65 


63.2 


3.71 


3.1 M 
1.6 ~J 




G038.94-00.46 


38.94 


-0.46 


41.6 


2.97 


4.33 


39.01 


-0.51 


0.36 


0.28 


-10 


0.23 


1.90 


31.8 


1.02 




G023.44-00.21 


23.44 


-0.21 


103.7 


3.44 


4.23 


23.55 


-0.27 


0.35 


0.23 


28 


0.12 


2.13 


51.7 


0.49 


2.5 




G030.79-00.06 


30.79 


-0.06 


94.7 


6.12 


4.23 


30.86 


-0.04 


0.39 


0.28 


66 


0.32 


2.24 


79.5 


2.68 


3.9 




G030.29-00.21 


30.29 


-0.21 


104.5 


3.01 


3.92 


30.36 


-0.18 


0.39 


0.23 


66 


0.26 


1.66 


40.0 


0.98 


2.0 




G053. 14+00.04 


53.14 


0.04 


22.0 


2.39 


3.88 


53.15 


0.09 


0.42 


0.26 


74 


0.29 


2.01 


26.8 


0.93 


1.3 




G022.44+00.34 


22.44 


0.34 


84.5 


2.81 


3.62 


22.51 


0.30 


0.41 


0.23 


32 


0.14 


1.52 


29.0 


0.42 


1.4 




G024.49+00.49 


24.49 


0.49 


102.4 


5.24 


3.58 


24.48 


0.30 


0.52 


0.32 


8 


0.46 


1.46 


67.8 


2.21 


3.3 




G049.39-00.26 


49.39 


-0.26 


50.9 


3.54 


3.54 


49.44 


-0.23 


0.30 


0.18 


63 


0.15 


1.93 


54.1 


0.90 


2.7 




G019. 39-00.01 


19.39 


-0.01 


26.7 


3.88 


3.48 


19.53 


0.02 


0.42 


0.26 


-2 


0.31 


1.73 


39.0 


2.00 


1.9 




G034.74-00.66 


34.74 


-0.66 


46.7 


4.33 


3.46 


34.88 


-0.63 


0.42 


0.34 


51 


0.34 


1.91 


35.3 


2.27 


1.7 


Y 


G023.04-00.41 


23.04 


-0.41 


74.3 


4.20 


3.40 


23.07 


-0.45 


0.50 


0.31 


7 


0.45 


1.54 


45.9 


2.31 


2.3 




G018. 69-00.06 


18.69 


-0.06 


45.4 


3.86 


3.33 


18.77 


-0.13 


0.26 


0.18 


-22 


0.13 


1.87 


32.8 


0.85 


1.6 




G018. 19-00.31 


18.19 


-0.31 


50.1 


4.15 


3.29 


18.20 


-0.44 


0.39 


0.27 


89 


0.27 


1.89 


52.5 


1.52 


2.6 


X 


G025.64-00.ll 


25.64 


-0.11 


93.9 


2.78 


3.21 


25.67 


-0.26 


0.36 


0.30 


12 


0.29 


1.64 


24.8 


1.06 


1.2 


Y 


G024.79+00.09 


24.79 


0.09 


110.4 


3.25 


3.17 


24.62 


0.17 


0.48 


0.37 


4 


0.48 


1.50 


50.3 


1.79 


2.5 





Table 3: Properties of the molecular clumps identified in the GRS. 



Cloud 


Clump 


Peak 


Vls.r 


A£ 


Ab 


AV 




A 


Ipeafc 


^-total 


N(H 2 ) 


Flag 


GRSMC 




t 


b 


























(°) 


(°) 


(kms -1 ) 


(°) 


(°) 


(kms -1 ) 


(K) 


(dcg 2 ) 


(Kkms- 1 ) 


(Kkms -1 dcg 2 ) 


(xlO 21 cm" 2 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


G053.59+00.04 


cl 


53.57 


0.06 


23.5 


0.08 


0.08 


1.95 


21.52 


2.00E-02 


35.65 


2.39E-01 


10.6 


V 


G053.59+00.04 


c2 


53.62 


0.03 


23.5 


0.21 


0.09 


1.49 


20.23 


4.55E-02 


34.02 


2.93E-01 


9.9 


- 


G053.59+00.04 


c3 


53.40 


0.03 


22.5 


0.14 


0.08 


2.39 


7.60 


1.34E-02 


9.82 


4.43E-02 


3.7 


X 


G053.59+00.04 


c4 


53.49 


-0.01 


21.6 


0.04 


0.07 


1.52 


5.98 


4.84E-03 


3.91 


8.80E-03 


2.9 


V 


G029.89-00.06 


cl 


29.87 


-0.04 


100.9 


0.07 


0.09 


2.52 


16.06 


2.20E-02 


46.52 


3.41E-01 


7.9 


- 


G029.89-00.06 


c2 


29.96 


-0.02 


97.7 


0.06 


0.10 


2.20 


15.75 


1.90E-02 


47.69 


1.91E-01 


7.7 


- 


G029.89-00.06 


c3 


29.91 


-0.04 


98.3 


0.08 


0.07 


2.43 


14.52 


1.73E-02 


46.81 


2.02E-01 


7.1 


- 


G029.89-00.06 


c4 


29.91 


-0.05 


101.5 


0.08 


0.10 


2.71 


14.21 


2.61E-02 


44.25 


3.52E-01 


7.0 


- 


G029.89-00.06 


c5 


29.93 


-0.02 


98.1 


0.04 


0.04 


2.95 


13.06 


5.90E-03 


36.76 


6.97E-02 


6.4 




G029.89-00.06 


c6 


30.01 


-0.04 


93.9 


0.05 


0.06 


2.80 


12.06 


1.30E-02 


27.99 


1.24E-01 


5.9 


V 


G029.89-00.06 


c7 


30.00 


0.01 


100.2 


0.08 


0.08 


3.04 


10.27 


2.07E-02 


26.91 


1.58E-01 


5.0 




G029.89-00.06 


c8 


29.78 


-0.26 


99.8 


0.11 


0.08 


2.02 


9.42 


2.24E-02 


19.33 


1.25E-01 


4.6 




G029.89-00.06 


c9 


29.69 


-0.17 


100.2 


0.10 


0.10 


2.32 


8.38 


1.98E-02 


19.14 


1.45E-01 


4.1 


X 


G029.89-00.06 


clO 


30.03 


-0.51 


100.5 


0.13 


0.06 


0.81 


7.96 


9.87E-03 


4.14 


1.61E-02 


3.9 




G029.89-00.06 


ell 


30.10 


-0.26 


96.0 


0.07 


0.09 


2.02 


7.44 


3.25E-03 


9.42 


1.11E-02 


3.7 




G029.89-00.06 


cl2 


30.09 


-0.29 


96.8 


0.03 


0.05 


1.62 


7.38 


2.80E-03 


7.50 


7.72E-03 


3.6 




G029.89-00.06 


cl3 


30.01 


-0.26 


101.7 


0.12 


0.17 


2.45 


7.21 


2.68E-02 


13.08 


9.75E-02 


3.5 




G029.89-00.06 


cl4 


29.69 


-0.22 


95.4 


0.07 


0.06 


2.79 


7.21 


2.31E-03 


9.57 


6.18E-03 


3.5 


X 


G029.89-00.06 


cl5 


30.09 


0.13 


97.1 


0.10 


0.08 


1.75 


6.83 


9.12E-03 


7.45 


2.38E-02 


3.4 


Y 


G029.89-00.06 


cl6 


30.31 


-0.28 


99.4 


0.06 


0.13 


0.55 


6.81 


2.80E-03 


4.94 


6.23E-03 


3.3 




G029.89-00.06 


cl7 


30.26 


-0.29 


100.9 


0.06 


0.04 


0.81 


6.65 


2.57E-03 


4.88 


4.70E-03 


3.3 




G029.89-00.06 


cl8 


30.04 


0.01 


104.5 


0.09 


0.08 


2.59 


6.31 


9.38E-03 


4.68 


1.22E-02 


3.1 


V 


G029.89-00.06 


cl9 


30.04 


0.04 


105.1 


0.11 


0.05 


2.00 


6.29 


2.12E-03 


3.10 


2.74E-03 


3.1 


V 


G029. 89-00.06 


c20 


30.13 


-0.33 


99.8 


0.09 


0.10 


1.40 


5.62 


6.24E-03 


6.51 


1.21E-02 


2.8 





Table 4: Summary of the derived properties of the GRS molecular clouds and clumps. 



Properties Clouds Clumps 





Min 


Max 


Mean 


Median 


Std dev 


Slope 


K-S test 


Min 


Max 


Mean 


Median 


Std dev 


Slope 


Peak T m (, temperature (K) 


0.8 


5.8 


1.6 


1.4 


0.7 


-3.1 ± 0.2 


0.000000 


1.2 


27.3 


5.2 


4.7 


2.0 


-4.0 ± 0.3 


Lincwidth (kms -1 ) 


0.6 


9.8 


3.6 


3.4 


1.4 


-2.9 ± 0.4 


0.001645 


0.2 


8.3 


1.4 


1.2 


0.8 


-3.1 ± 0.4 


Scmimajor Axis (dcgs) 


0.06 


1.16 


0.41 


0.41 


0.15 


-5.3 ± 0.3 


0.000014 


0.01 


0.25 


0.03 


0.02 


0.02 


-2.3 ± 0.2 


Semiminor Axis (dcgs) 


0.05 


0.53 


0.23 


0.22 


0.09 


-2.7 ± 0.5 


0.000138 


0.01 


0.24 


0.03 


0.02 


0.02 


-2.0 ±J.3 


Position Angle (degs) 


-90.0 


90.0 


1.2 


-1.0 


43.5 




0.965160 














Log[Peak H2 column density (cm -2 )] 


20.7 


23.0 


21.8 


21.8 


0.3 


-2.3 ± 0.2 


0.000000 


20.8 


22.1 


21.4 


21.4 


0.2 


-4.1 ±10.2 


Excitation Temperature (K) 


4.1 


17.5 


8.8 


8.5 


2.0 


-6.8 ± 0.6 


0.000001 














Opacity 


0.07 


0.50 


0.13 


0.12 


0.04 


-3.7 ± 0.3 


0.000231 














Radius (pc) 


1.6 


97.5 


24.1 


22.6 


12.5 


-4.0 ± 0.4 


0.002011 


1.2 


27.3 


5.2 


4.7 


2.0 


-4.1 ± 0.2 


Log[LTE Mass (M©)] 


2.2 


5.7 


4.5 


4.6 


0.6 


-1.8 ± 0.3 


0.000000 


0.2 


5.3 


2.9 


2.9 


0.7 


-1.2 ± 0.1 


Log[Virial Mass (M©)] 


1.9 


5.7 


4.5 


4.5 


0.6 


-1.5 ± 0.2 


0.000532 


0.3 


5.2 


2.7 


2.7 


0.7 


-1.2 ± 0.1 


Log[Density (cm -3 )] 


0.6 


2.8 


1.7 


1.7 


0.4 


-1.4 ± 0.2 


0.000000 


0.8 


3.6 


2.5 


2.6 


0.4 


-2.9 ± 0.4 



